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1 Introduction 

The theory of derivatives of non-integer order goes back to Leibniz, Liouville, Riemann, Grun- 
wald, and Letnikov P3I21E1- Fractional analysis proved to be useful in mechanics and physics. 
In a fairly short time the list of such applications continuously grows. The applications include 
chaotic dynamics jHEl, material sciences jHllZllHlE], mechanics of fractal and complex media 
[inilIIllI2], quantum mechanics [131 Ej, physical kinetics jHHSllinilHl, plasma physics [111111], 
electromagnetic theory [201 1^ long-range dissipation |22], non-Hamiltonian mechanics 
[231 121] , long-range interaction [23 [20] , anomalous diffusion and transport theory [H |2Z1 12H1 121] • 

In this paper, we consider the fractional generalization of nonholonomic constraints such 
that the constraint equations consist of fractional derivatives, called fractional constraints. The 
corresponding equations of motion will be derived by the d'Alembert-Lagrange principle and 
some simple examples are considered. 

In Sec 2, we provide a brief review of nonholonomic systems, fix notations and convenient 
references. In Sec. 3, we consider fractional generalizations of nonholonomic constraints. Some 
examples are considered. In Sec. 4, we discuss the applicability of the stationary action principle 
for fractional constraints. Finally, a short conclusion is given in Sec. 5. 
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2 Nonholonomic Constraints with Integer Derivatives 



In this section, a brief review of nonholonomic systems is considered to fix notations and provide 
convenient references j32j . 



2.1 Lagrange Equations for Nonholonomic System 

It is known that the d'Alembert-Lagrange principle allows us to derive equations of motion 
with holonomic and nonholonomic constraints. For N-particle system it has the form of the 
variation equation 

^ -F, Ur. = 0, 



dt 

where {i = 1...N) is a radius- vector of ith particle, Vj = !■« is a velocity, and Fj is a force that 
acts on ith particle, and the sum over repeated index i is from 1 to A^. To exclude holonomic 
constraints, the general coordinates {k = 1, ...n) are used. Here, n = 3N — m is a number 
of degrees of freedom, where m is a number of holonomic constraints. Then is a function of 
generalized coordinates and time: = ri{q,t). Using 5rj = {dri/dq^)5q^ , Eq. (P) gives 

and sum over repeated index k is from 1 to n. Then, we define j^Dl the generalized forces: 

dv 

Qfc = Fi7r4 A; = l,...,n. 

oq^ 

By usual transformations jSOl 

dimwi) dvi d f dri \ d dvi 

- — rriYi—r - {m\i' - 



dt dq'' dt \ dq'' J dt dq 



d f dwi\ dvi 
— Tn\i -—- — mv,- 



dt \ dq'' J dq'' 



d d i"m \ d /m \ d dT dT 

ViV,- - —r — ViVi 



dt dq>'\2 ' 7 dq''\2 ' 7 dt dq'' dq'' 
we transform Eq. Q into 



^ dt dq'' dq'' 
where T = is a kinetic energy. Using 

dvi{q, t) dvi dq'' dr^ 



d dT dT ^ ^ , k n 

- T7T - Qfc = 0, (3) 



dt dq'' dt dt ' 



+ ^> (4) 
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we get 

T = - {gki{q, t)q'q' + 2g,{q, t)q' + g{q, t)) , 

where 

dridvi dvidvi dvidvi 

9Uq,t) = ^^, 9u{q,t) = —^-, 9iq,t) = (5) 

For the nonholonomic constraint, 

RkSq" = 0, (6) 
where Rk is a reaction force of the constraint 

/(g,g) = o, (7) 

and the variations Sq'^ are defined [HIl IS2] by 

P^'-O. (8) 
Equation (jSI) is called Chetaev's condition [32] • Comparing Eqs. and (jS)), we obtain 



df 

= (9) 



where A is a Lagrange multiplier. Chetaev's definition of variations states that the actual 
constrained motion should occur along a trajectory obtained by normal projection of a force onto 
a constraint hypersurface. The constraint force Rk is minimum when is chosen perpendicular 
to the constraint surface or parallel to the gradient df /dqk- 

In general, the nonholonomic system is subjected to action of the generalized force Qk and 
the constraint force Rk- Then the variational equation is 



Qk - Att^ Sq' = 0. (li: 



From 0, we obtain 

dt dq^ dq^ ^ dq^ _ 

In Eq. (jSJ, we can consider 5(f, s = 1,2, ...,n — 1, as independent variations. Then 5q^ is not 
independent, and Eq. (jHl) gives 



' s=l 



Suppose that A satisfies the equation 



d dT dT df 

JtW'W^'^'^W^^' ^^^^ 
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Then the term with k = n in (fTT|) is equal to zero, and Eq. (fTT|) has n — 1 terms with 
k = 1, ...,n — 1. In Eq. (jllll . the variations with k = 1,2, ...,n — 1 are independent, and the 
sum is separated on n — 1 equations. As the resuh, Eq. (fTTj) is equivalent to 

d dT dT df 

dtW-W = ^'^^W^ (^ = l,...,n). (13) 

Equations ((Tj) and ()13|) form a system of n + 1 equations with n + 1 unknowns A and q'^ , where 
k = 1, ...,n. Solutions of these equations describe particles motion as a motion of system with 
nonlinear nonholonomic constraint ((Tj). 

2.2 Nonholonomic System as a Holonomic One 

In this section, we present the equations of motion with nonholonomic constraint as equations 
for a holonomic system. 

The canonical momenta p'^ are defined by 

dT 

Pk = -gT^ = mgkiiq,t)q^ + mgkiq,t), {k=l,...,n). (14) 
Using (^H), we can define 

f{p,q,t) = fiqiq,p.t),q,t). (15) 

Suppose that the constraint is integral of motion. Then the total time derivative of ()15|1 gives 

-^ = +^•'^ + ^ = (16) 

dt ' dpk dq'^ dt 

Substitution of ^ into ^ get 

^i— + Q, + \^ +^q' + ^ = 0. (17) 
dpk \ dq^ dq^ j dq^ dt 

From Eq. (fTTj) . we obtain 

df dfY' ( df fdT ^ ^ \ ^ df , df 



dpk dq^ ] \ dpk \dq'' ~'~ ^'^J ~'~ ~'~ dt 



Then the Lagrange equations ()T3|l have the form 



-1 



d dT dT ^ df I df df \ ( df f dT ^ ^ \ ^ df ^ df 



dtdq^ dq^ dq^ ydpkdq'^j ydpk \(9g^ + '^fc J + q^^Q + 1 • (^9) 

Equations (fT^ describes the motion of a holonomic system with n degrees of freedom. For any 
trajectory of the system in the phase space, we have / = 0. If the initial values gfc(O) and gfc(O) 
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satisfy the constraint condition /(^'(O), g(0), to) = 0, then the solution of Eq. ^T^i is a motion 
of the nonholonomic system. 

Let us define a generahzed force Ak = Qk + Rk, which depends on generahzed velocities q'^, 
generalized coordinates q^, and time t. If 

dAk ^ dAm ^ Q 

Qqm Qqk ~ ' 

dAk dKn^ld_ ^dAk _ dKA 
Qqm Qqk 2dt ydq"^ dq^ ) ' 

known as the Helmholtz conditions, are satisfied, then a generalized potential U = U{q, q, t) 
exists and 

d dU dU _ 
dt dq'' dq'^ ^ 

In this case, the Hamilton variational principle has the form of the stationary action principle. 
In order to use this principle for a nonholonomic system, we should consider such trajectories 
that their initial conditions satisfy the constraint equation ((7j). 

Note that nonholonomic constraint ((7j) and non-potential generalized force Qk can be com- 
pensated such that resulting generalized force is a generalized potential force, and system 
is a Lagrangian and non-dissipative system with holonomic constraints. 



3 Constraints with Fractional Derivatives 
3.1 Fractional Derivatives 

The fractional derivative has different definitions [Ij, |2j|, and exploiting any of them depends 
on the kind of the problems, initial (boundary) conditions, and the specifics of the considered 
physical processes. The classical definition is the so-called Riemann-Liouville derivative [H |2] 

.mi.) ' ^ r m"^ 



T{m-a) dx"^ (x - 2)"-™+i ' 

where m — 1 < a < m. Due to reasons, concerning the initial conditions, it is more convenient 
to use the Caputo fractional derivatives [71 ESI EE] • Its main advantage is that the initial con- 
ditions take the same form as for integer-order differential equations. 



Definition. The Caputo fractional derivatives are defined by the equations 

r(m -a) {t-T 



1 /■* f('"Hr) 
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''■'^ ' r(m - a) Jt (t - r)°+i-™ ' ^ 



where m — 1 < a < m, and /'•'"■' (t) = d^f{T)/d 



T" 



Proposition 1. The total time derivative of the Caputo fractional derivative of order 
can be presented as a fractional derivative of order a + 1 by 

J j-m—a—l 

where m = [a] + 1, and [...] means floor function. 

Proof. The definition (j?H) can be presented in the form 

aDtf = air^'DTf, (2^ 

where aJr~'^ is a fractional integral 



The operations and do not commute: 

Dl aJtfit) =a JtD]f{t) + ^/(a). (2C 

From (j24p . we have 

I aDtf{t) = D] aJlDTf{t) = D] aJ^f^-'Kt), (2< 

where aJf is a fractional integration of order e = m — a. Using and (jTTj) . we get 

D]JtDTf{t) =a JtD]f^^\t) + ^/^'"Ha) = 

Substitution of (PHj) into ^F!\ proves 



3.2 Fractional Equations of Motion 

Asume that the constraint equation has fractional derivatives: 

f{q,q, aDtq, tD^q)=0, (29) 

i.e. it is a fractional differential equation 0. Such constraint can be called fractional nonholo- 
nomic constraint. Since Eq. ()29p has also derivatives of integer order, we can use the Chetaev 
definition of variation (jHJ and the Lagrange equations ()13|) . For generalized potential forces 

ddU dU 

^''^ ^f7^ ~ 7^' (fc=l,...,n), 
at oqk oqk 

and we can rewrite Eq. (fTSj) as 

d dL _d£^ ^9/ ^ ^ ^ 
dtdqk dqk dq^ 

where L = T{q, q) — U{q, q) is the Lagrangian. To simphfy our calculations, we consider 

L = L{q,q) = ^-j2i4kr-u{q), (31) 

k=l 

where u{q) is a potential energy of the system. Then, Eq. (jHUj) becomes 

du f 

Suppose that the constraint is an integral of motion, i.e., df/dt = 0. Then 

dldqk df djaPM df djtPM dldqk ^ 

dqk dt d{aDfqk) dt ditD^qk) dt dqk dt ' ^ ' 

Equation (jH^ can be presented as 

IT^^ + % ^ aDtqk + M ,D^qk + = 0. (34) 

dqk d[aDfqk) d{tD^qk) dqk 



Substitution of ^ into (j3D) gives 

df f du ^ ^ df \ ^ df 2 

, 9f , df 



I +^7^)+ .^f na ^ D] aDtqk + ^, ' M ,D^qk + -^qk = 0- (35) 

dqk V 9qk dqk J d[aDfqk) d{tD°qk) dqk 

From this equation, one can obtain the Lagrange multiplier A: 

df df f9f du df 



X 



dq'm dqm ) V 9<li 9{aDfqi ] 



D] aD^qi - ^, y , D] tD^qi - ^qA . (36) 
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Insertion of Eq. (jHHjl into Eq. (j!?^ yields 



dq,^ dq,{dqrr.dq^) [dq, dq, ^UA^g,) * ' ^GD.^gz) * * ^g.^J " 

(37) 

These equations describe holonomic system that is equivalent to the nonholonomic one with 
fractional constraint. For any motion of the system, we have / = 0. If the initial values satisfy 
the constraint condition f{q{0),q{0), aD^q{0), tD^g(O)) = 0, then the solution of Eq. (jHTj) 
describes a motion of the system (jHT|) with fractional constraint 



3.3 Linear Fractional Constraint 

Suppose that the constraint (pUj) is linear with respect to integer derivatives g^, i.e., 

f = akqk + f3{aD'^q, tD^q,q). (38) 
In this case, Rk = a^, and Eqs. (|H7|) can be presented as 

-A / akai\ du flfc-A/ df ^ df ni n" , • ^ 

1k = -y[Oki T- la 2 na qi + tDf^ qi + —qi , 

^ V a2 / dqi a2 \d{aDfqi) d{tD^qi) dqi J 

(39) 

where = J2k=i ^kO'k- If 

/?(,A"g, tDtq,q)=huaDtqk. (40) 

then 

/ = flfegfc + hk aDtqk. (41) 

This constraint is linear with respect to integer derivative and fractional derivatives aD^qk- 
Then the equations of motion are 

= - E [6k, - ^ j ^ - E .A ^^ (42) 

1=1 ^' i=i 

Using Proposition 1, we obtain 

where q^"^\a) = {D'^q{t))t=a- As the result, we get the fractional equations of motion with 
Caputo derivative of order a + 1. The nonholonomic systems with integer a are considered 
inRefs. ^Mli- 
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3.4 One-dimensional Example 

In one-dimensional case {n = 1), Eq. (jl^ with oD^q has the form 

q = -^Dl ,Dtq. (44) 

Eq. can be presented as 

A'[g+(&iM)oA"g]=o. (45) 

Then 

g + (6i/ai)oA"? = Co. (46) 
Supposing a > 1, and using proposition 1, we get 

Dl[q + {h/a,)oDr'q] = — Tg(0) + Co- (47) 

Oii [m — a + I) 

As the result, we obtain 

j.m—a+1 

oDr'q + {ai/h)q = — -g(0) + C,t + C^, (48) 

L [m — a + 2) 

where we use xT{x) = T{x + 1), and Ci = Coai/bi. 

For 2 < a < 3, Eq. PH|) describes the linear fractional oscillator 



oDr'q{t)+u;^q{t) = Q{t), (49) 
where uj'^ = (ai/bi) is dimensionless "frequency", and Q{t) is the external force: 



j.m— o+l 



T{m — a + 2) 

The Caputo fractional derivative o-D"~^ allows us to use the regular initial conditions jS] for 
Eq. (jini). The linear fractional oscillator is an object of numerous investigations \?)7 \ 138 } l39 | 
mi lU ia ESI in iSl EE! because of different applications. 
The exact solution |SZl EB] of Eq. (gHI) for 2 < a < 3 is 

g(t) = g(0)E,_i,i(-^V-i)+tg'(0)E„_i,2(-^V-^)- / Q{t - r)qo{r)dr, (50) 

^0 



where 



OO J, 



is the generalized two-parameter Mittag-LefHer function 53EB|; and 

go(r)=^,_i,i(-^V-i). 
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The decomposition of (|Hn|l is ;37i| : 

q{t) = g(0) [fa,o{t) + 9a,o{t)] + tq{0) [U,i{t) + g^,i{t)] - [ Q{t - T)qo{T)dT, (52) 

where 







^ Jo + 2r'^cos(7ra) + 1 

ga kit) = ^e*™'("/") cos [t sin(7r/a) - nk/a] , [k = 0, 1). (53) 
a 

For the initial conditions g(0) = 1, and g(0) = 0: 

g(t) = K(-t") = /«,o(t) + Qafiit) - [ Qit- r)[/,,o(r) + g.fi{r)]dT. (54) 
The first term in ()54p decay in power law with time while the second term decays exponentially 

EH EHi Ennui. 

3.5 Two-dimensional Examples 

In the two-dimensional case (n = 2), Eq. ()42|) has the form 

al du aia2 du Qi&i ^i^ni . 

«i du aiaa .9n ^2^1 1 „ ^2^2 1 „ . . 

Let us consider the special cases of these equations. 



(1) Suppose Oi = 0; then 



qi = q2 = --D] aDUi - -D] aD^q^. (57) 

dqi 02 (22 



If ai = 0, and 62 = 0, then (|K7jl are 



c'gi 02 



(2) Suppose 61 = and ai = 02 = c; then we have 

1 du 1 du 62 



^^ = -2 5^ + 2 5^-2^^-^^^^- 
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Using 

we can rewrite Eqs. (I59p and in the form 



gi + ^2 qi-q2 , , 

X= , y= , 9 = 02/0, 



dU 

x = -gDl aDtx + gD] aDty, Z/ = (61) 

where U{x, y) = u{qi, ^2) = u(yX + y,x — y). If f/ = K{x)y + s{x), then Eq. is 

X = -gD] aD^x + gD] aD^y, y = ~K{x). (62) 
Using Dfy = Df'^y, Eq. ^ gives 

X = -gD] ^Dtx - gD] aDr^K{x). (63) 

Then 

Dl[x + g aDtx + g aD^^K^x)] = 0. 

For a > 2, then we can get 

X + g-^D'^x + gaDr^K{x) + C = 0, (64) 

where C is a constant. Using aJtfif) = /(^)) ^"^^ aD^ aD^f(t) = D^f{t), where e = m — a, 
Eq. ()64p can be written as 

D'x + ^-^x^™) + gaD'^-^K{x) = 0. (65) 

If m = 2 (1 < a < 2), then 

g-^i + D'x + gaK{x) = 0. (66) 

This equation can be considered as an equation of nonlinear fractional oscillator jHH HT] , where 
fractional derivative describes the power dumping. 

4 Fractional Conditional Extremum 

4.1 Extremum for Fractional Constraint 

Let us consider the stationary value of an action integral 



6 f dt L{q, q) = 0, 

J a 



for the lines that satisfy the constraint equation f{q,q) = 0. Using the Lagrange multiplier 
/X = /i(t), we get a variational equation 



5 I dt[Liq,q)+ixfiq,q)] = 0. 

J a 
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Then the Euler-Lagrange equations are 

ddL dL (dS ddf\ .9/ ,. , , 

Note that these equations consist of the derivative of Lagrange multipher /i. The proof of Eq. 
(I67|l is reahzed in [34]. 

For the fractional constraint 

aDU. tD^q)=0, (68) 

we can define the Lagrangian as 

Ciq, q, aD^q, tD^q, A) = L{q, q) + g, ^D'^q). (69) 

Using the Agrawal variational equation [HBl, we obtain the Euler-Lagrange equations 

^ _ + - (k-l n) (70) 

dq, dtdq^^^^'dUDM^'^'di^DM' ' ik-l,...,n). (70) 

Substitution of Eq. ^ into Eq. ^ gives 

dL d dL Of ^^f Of \ ^„ / df \ d f df\ ^ 



dqk dtdqu ' y^d{aD'^qu)J " ditD^ J c^H -^g, 

These equations describe the fractional conditional extremum. 

Let us consider applicability of stationary action principle for mechanical systems with frac- 
tional nonholonomic constraints. The equations of motion are derived from the d'Alembert- 
Lagrange principle. The fractional conditional extremum can be obtained from stationary 
action principle. In general, these equations are not equivalent [21] • The condition of this 
equivalence for fractional constraints is suggested in the proposition. 

Proposition 2. Equations A3U^} and ^71\ ) for nonholonomic system with fractional con- 
straint 116^} have the equivalent set of solutions if the conditions 



are satisfied. 



Of 

5qk = 0, -^5gfc = 0. (72) 
dqk 



Proof. To prove the proposition, we multiply Eqs. (jHUj) and ()71|) on the variation 5q^ and 
consider a sum with respect to k: 

ddL dL\ df 

at oqk oqk / oqk 
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r dL _ d dL\^ ^ 
\dqk dtdqkj 



' df ( _df_\ ( df \ _d ( df 



From the definition of variations (jS)), Eq. (f73|) is 

d dL dL\ 



Sqk = 0, 
(74) 



dtdqt dqkj 
Substituting Eq. ^ into Eq. (HH), we obtain (17^ . 



6qk = 0, (75) 



It is known tliat stationary action principle cannot be derived from the d'Alembert- 
Lagrange principle for wide class of nonholonomic and non-Hamiltonian systems. The same 
can be applied to the case of nonlinear fractional nonholonomic constraints. 

4.2 Hamilton's Approach 

Using Eq. (f70|) . we define the momenta 

dC dL df 

Pk = TT^ = f^TT^, (76) 

dqk dqk dqk 

and the Hamiltonian 

n{q,p) = Pkq - C, (77) 

where C = L + fif . Eq. (fTOj) gives 

To simplify our calculations, we consider the Lagrangian 

L = l{qr-u{q), 

and the fractional nonholonomic constraint 

/ = Ak{q, aDtq) qu = 0, {a^l). (79) 
From Eqs. (|76|) and (|79|) . we obtain 

Pk = qk + ^iAk{q, aDtq). (80) 

Then the Hamilton equations are 

qk=Pk- fJ'Akiq, aDfq), (81) 



13 



du{q) , , r.af..,dAi{aD^q) 



oqk oqk \ OaDfqk J 

To find the Lagrange multiplier /i = /i(t), we multiply Eq. (jSUj) on the functions and consider 
the sum with respect to k: 

AkPk = Akqk + fJ'AkAk = (83) 

Here, we use the constraint and the notation A"^ = AkAk, where Ak = Ak{q, aDfq). From 
(JHSl), we get 

AkPk 

Substitution of ^ into Eqs. dHU) and (jHS) gives 

% = (^Ski - Ph (85) 

du{q) AmPm dAi [ AraPm dAi . \ 



If Ofc = 0, then we have usual equations of motion for Hamiltonian systems. Note that we derive 
Hamiltonian equations from Euler-Lagrange equations without using the Legendre transforma- 
tion, which is typically used. 

5 Conclusion 

The classical mechanics of nonholonomic systems has recently been employed to study a wide va- 
riety of problems in the molecular dynamics |49j. In molecular dynamics calculations, nonholo- 
nomic systems can be exploited to generate statistical ensembles as the canonical, isothermal- 
isobaric and isokinetic ensembles [SniinillS21ISSllSlEniEniEZllSHlEniEni- Using fractional 
nonholonomic constraints, we can consider a fractional extension of the statistical mechan- 
ics of conservative Hamiltonian systems to a much broader class of systems. Let us point out 
some nonholonomic systems that can be generalized by using the nonholonomic constraint with 
fractional derivatives. 

(1) In the papers [201 Ol E21 E3] j the constant temperature systems with minimal Gaussian 
constraint are considered. These systems are non-Hamiltonian ones and they are described by 
the non-potential forces that are proportional to the velocity, and the Gaussian nonholonomic 
constraint. Note that this constraint can be represented as an addition term to the non-potential 
force [UUj . 

(2) In the papers [HHlEni, the canonical distribution is considered as a stationary solution 
of the Liouville equation for a wide class of non-Hamiltonian system. This class is defined by 
a very simple condition: the power of the non-potential forces must be proportional to the 
velocity of the Gibbs phase (elementary phase volume) change. This condition defines the 
general constant temperature systems. Note that the condition is a nonholonomic constraint. 
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This constraint leads to the canonical distribution as a stationary solution of the Liouville 
equations. For the linear friction, we derived the constant temperature systems. A general 
form of the non-potential forces is derived in Ref. [UU] . 
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